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Abstract 

The internal structures and compositions of Uranus and Neptune are not well constrained due to the 
uncertainty in rotation period and flattening, as well as the relatively large error bars on the gravitational 
coefficients. While Uranus and Neptune are similar in mass and radius, they differ in other physical prop- 
erties such as thermal emission, obliquity, and inferred atmospheric enrichment. In this letter we consider 
the uncertainty in the planetary rotation periods, show that rotation periods more consistent with the 
measured oblateness imply that Uranus and Neptune have different internal structures, and speculate 
on the source of that difference. We conclude that Uranus and Neptune might have very different struc- 
tures and/or compositions despite their similar masses and radii. We point out that understanding these 
differences can have important implications for our view of the formation and evolution of Uranus and 
Neptune as well as intermediate-mass extra-solar planets in general. 

1 Introduction 

Uranus and Neptune are often thought of as twin planets. They both formed in the outer solar system, 
and their masses and radii are very similar, as are their gravitational moments. Even their rotation periods 
are within 10% of each other. Yet, there are important differences between these planets. Uranus' mass 
is slightly smaller than Neptune's, but its radius is somewhat larger, so the difference in mean density, is 
considerable. Indeed interior models of Uranus and Neptune (Podolak et al. 1995; Helled et al. 2011) 
indicate that Neptune's outermost envelope is more enriched in high-Z material than that of Uranus. 

In addition, Neptune has an internal heat source, while Uranus is in equilibrium with solar insolation 
(Pearl et al. 1990; Pearl & Conrath 1991) suggesting that Uranus' interior may not be fully convective, 
and/or that it contains compositional gradients which hinder convection. The difference in thermal flux is 
indicated from calculations of Uranus' thermal evolution which find that unlike Neptune, Uranus cannot 
reach its measured intrinsic luminosity within solar system's age if an adiabatic interior is assumed (Fortney 
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& Nettelmann 2010, Fortney et al., 2011). Another distinct feature of Uranus is its large axial tilt. This 
was almost certainly caused by a dramatic event in the early history of the solar system (Safronov 1966), 
although an alternative explanation for Uranus' tilt is orbital migration (Boue & Laskar, 2010). Stevenson 
(1986) suggested that this same dramatic event might have caused the necessary compositional gradients for 
inhibiting convection. We return to this point later. 

Despite the available measurements, the internal structures of Uranus and Neptune are not well con- 
strained, and as a result, their bulk compositions are essentially unknown (Podolak et al. 1995; Marley et 
al. 1995; Helled et al. 2011). In fact, it is still unclear what the mass fraction of water is in Uranus and 
Neptune, despite their categorization as 'icy planets'. 

Recently, Nettelmann et al. (2012), showed that by assuming modified solid-body rotation periods for 
Uranus and Neptune, the derived internal structures can change considerably. It was concluded that an 
uncertainty in the planetary rotation period can be crucial for constraining the internal structure. In this 
letter, we investigate further how the uncertainties in the gravity field, shape, and rotation period might 
affect the inferred planetary interior. 

2 Observed Parameters 

Interior models use observed physical properties to constrain the planetary composition and its depth de- 
pendence. Below we summarize the measurements and suggest that their uncertainties, when included in 
interior models, can be crucial for understanding the formation, evolution and interiors of the planets. 

2.1 Gravitational Moments 

The gravitational moments are probably the best known physical properties of Uranus and Neptune after 
their masses. Although the error bars are fairly large and only J 2 , J4 are known, the gravity data provide 
tight constraints on the planetary density profile. Recently, Jacobson (Jacobson 2007, 2009) rc-dctcrmined 
the gravitational harmonics Uranus and Neptune and provided more accurate estimates for the gravitational 
harmonics with significantly smaller error bars. The updated gravitational data of Uranus and Neptune are 
given in Table 1. 

2.2 Oblateness 

The planetary oblateness (flattening) is defined by 

Re-Rp 
~ Re 

where Re and Rp are the equatorial and polar radii, respectively. While knowledge of the planetary con- 
tinuous shape is available, typically, interior models use only the flattening. Although stellar and ring 
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occultations provide good determinations of the planetary shape (French et al. 1987), these data correspond 
to low pressure-levels and might not apply to the 1-bar level that is used by interior modelers since atmo- 
sphere dynamics (winds) can change the isobaric shape between pressure levels. The connection between 
the planetary shape at 1-bar and at lower pressure levels isn't entirely obvious and should be investigated in 
future research. In addition, there are no direct measurements of the polar radii for Uranus and Neptune, 
and they are inferred from the measured flattening at lower latitudes. 

Lindal ct al. (1987) and Lindal (1992) provided estimated the polar and equatorial radii of Uranus and 
Neptune using occultation data and these values were adopted as the "standard" for interior models (e.g., 
Guillot 2005). These shapes, however, were not inferred from the Voyager 2 radio periods but under the 
assumption that the planets rotate differentially on cylinders. This explains the inconsistency between the 
typically used flattening of Uranus and Neptune and the Voyager 2 rotation periods (Helled et al. 2010, 
hereafter HAS10). If differential rotation is confined to a thin atmospheric layer, and the deeper interior 
rotates uniformly, the shape of the geopotential surface is determined by this latter rotation period. For 
interior models the shape of the geoid should be used, since it represents most of the mass, but this shape 
is essentially unknown. 

2.3 Rotation Period 

The rotation periods of Uranus and Neptune are determined from radio and magnetic field data. HAS10 
suggested that, based on our experience with Saturn (Gurnett et al. 2007, 2010), the measured radio periods 
of Uranus and Neptune, 17.24 and 16.11 hr, respectively, might not represent the rotation periods of the 
deep interior. In addition to radio data, Voyager 2 also provided magnetic field data. For Uranus, Ness et al. 
(1986) obtained a rotation period of 17.29 ±0.1 hr; a period which is close to the radio period. For Neptune, 
the time base for the flyby was insufficient to determine the period accurately but a rotation period of <~ 16 
hr was found to fit the magnetic field data fairly well (Ness et al., 1989). Recently, Karkoschka (2011) derived 
a rotation period of 15.9663±0.0002 for Neptune by following stationary atmospheric features. This period 
is only slightly shorter than the Voyager radio period. However, the rotation period of the atmosphere might 
not be identical to the rotation of the deep interior and therefore should be taken with caution. Finally, 
the complex nature of the magnetic fields of Uranus and Neptune, and the fact that the location at which 
the magnetic field is generated, which could be at ~ 0.7 R p i an et (Stanley & Bloxham 2006), suggest that the 
magnetic period might not represent the rotation period of the bulk interior. 

Although the formal error bars on the measurements of angular speed, w, and / are fairly small, obser- 
vational ambiguities make these numbers less certain than the error bars indicate. These two parameters 
are used together with J 2n , for inferring the internal structures of Uranus and Neptune, and through them 
the structures of ice giants among the newly discovered exoplanets. It is therefore desirable to investigate 
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the implications of possible changes in these parameters since this may have important consequences for our 
understanding, not only of Uranus and Neptune, but also of a whole class of exoplanets. 



3 The Radau Approximation 



Radau (1885) introduced a useful approximation: if e = y/ f(2 — /) is the ellipticity of an equipotential 
surface in the rotating body, and a is the surface's mean radius, the Radau parameter is 

a de 



V = 



da 



Using this parameter and Clairaut's equation, one can show that 

d 



da 



pa b y/l + 1] =5pa 4 tp(r)) 



where, if p(a) is the density on an equipotential surface, 

ps if pr2dr 

is the average density of the enclosed volume, and 



The Radau approximation assumes ip(rj) = 1. 



3.1 Relations Between Observed Parameters 

The Radau approximation relates the moment of inertia to other, directly measurable, parameters. If M 
is the planet's mass, a is its mean radius, and G is Newton's constant, we can define a rotation parameter 
q = uj 2 R 3 /GM which is the ratio of the centrifugal to gravitational force. Using the Radau approximation 
the moment of inertia, C is related to q and / via (Zharkov & Trubitsyn 1978). 
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For a more precise version of the approximation (Yodcr 1995) we set g(r) = GM(r)/r 2 and define the 
constants 5\ and 82 by 
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The relation between the inertia factor and the observable parameters is given by 
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Note that when the Radau approximation is exact (ip(rj) = 1) S\ = 0. This isn't precisely true for the 
solar-system giant planets, but for simplicity we assume S\ = 0. 

The combination of gravitational and rotational forces causes an equipotential surface to be non-spherical. 
The requirement that the potential on a level surface be independent of latitude relates J2, /, and q. Yoder 
(1995) gives 



J4 ~ 



15 2 5 
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(5) 
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Eq.Q uses the flattening, which has traditionally been regarded as the least reliable of the measured pa- 
rameters. The usual practice among interior modelers is to assume a solid-body rotation period and use it 
to match the models to J2 and J4. However Eqs. (5][7) can be used to replace either q or / by a combination 
of J 2 and J4. 

The accuracy of the Radau approximation can be gauged from Fig.[l] Plotted is the inertia factor vs. 
rotation period as determined from one pair of observed parameters: q and / (blue), q and Ji (maroon) and 
J 2 and / (green). The solid and dotted curves are for Jupiter and Saturn, respectively. The observational 
parameters used are listed in Table 1, while the inertia factors computed for the different planets and rotation 
speeds are given in Table 2. 

Since the observational parameters for Jupiter and Saturn are well known, we would expect the three 
estimates for the inertia factor to coincide at the correct rotation period (shown by the vertical black lines) . 
As can be seen from the upper panel in the figure, the Radau approximation works extremely well for these 
two planets. For Jupiter it's only 35 s shorter than the system III period. For Saturn it's 7.72 m shorter 
than the Voyager radio period, but only 54 s shorter than the period determined by Anderson & Schubert 
(2007). It also suggests that Saturn is more centrally condensed than Jupiter which is consistent with detailed 
interior models (e.g. Guillot 2005). However, when the same quantities are plotted for Uranus and Neptune 
the disagreement between the period inferred from the Radau approximation and the Voyager radio period 
is striking as can be seen from the middle panel. Indeed HAS10 pointed out that the Voyager periods may 
be in error, and suggested rotation periods of 16.57 and 17.46 hr for Uranus and Neptune, respectively. 

When we use the Voyager radio periods for q, and fit Ji then, as can be seen from the middle panel in 
Fig.[T] Uranus and Neptune are found to have inertia factors 0.22 and 0.23, respectively, and their interiors 
are expected to be rather similar. On the other hand, when we use HAS10 rotation periods the agreement 
is much better. This is shown in the bottom panel of Fig. [l] which is similar to the bottom panel but with 
HAS10 rotation periods. The three approximations cross at 16.11 hr for Uranus and 17.17hr for Neptune. 
HAS10 rotation periods are not precisely equal to those determined from the Radau approximation but 
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clearly in better agreement. With the modified rotation periods the inertia factor for Uranus is ~ 0.22, while 
the value for Neptune is ~ 0.25 implying that Uranus is more centrally condensed than Neptune. Even the 
Voyager periods implied a slightly more condensed Uranus, in the sense that it's envelope was closer to solar 
composition, but the HAS10 periods greatly enhance this difference. 

The conclusion that Uranus is more centrally condensed can also be derived from the gravitational 
moments themselves. The higher order moments are more sensitive to the density distribution in the outer 
layers of the planet. For planets like Jupiter, where the core is a relatively small fraction of the volume, its 
contribution to both Ji and J4 is negligible, and the core properties are determined mostly from Jo (i.e., the 
mass). 

For Uranus and Neptune the central mass concentration occupies a large enough volume to contribute 
directly to J2 (Fig. (1) in Helled et al. 2011). In some models it is large enough to contribute to J4 as well. 
As a result, the ratio | J4/J2I contains information on the core's size. We ran models of Uranus and Neptune 
assuming that they consist of a constant density core surrounded by an envelope where the density is given 
by a 6th order polynomial in radius (e.g. Schubert et al. 2012). The core density is set to 10gcm~ 3 , which 
is appropriate for rocky material under the expected conditions. The polynomial coefficients are chosen so 
that the radius, mass, and J2 are reproduced to within the measured uncertainties. J4 is then computed 
for the model. The results are shown in Fig. ([2]), where the ratio 7 = — J4/J2 is plotted vs. the core radius 
(normalized to the planet's radius). 

These models show that 7 increases with increasing core size. For the Voyager periods (solid) Uranus' 
core can be as large as ~ 13% of the planetary radius. Neptune's core is less constrained because of the larger 
error bars on the measured moments, and may be substantially larger. For HAS10 periods (dotted) Uranus' 
core size is essentially unaffected, while Neptune's core must be even larger. In fact, for these models, there 
is actually a lower limit to the size of Neptune's core. While detailed conclusions about the core size cannot 
be derived from these simple models, they do show how the interior structure is affected. Recent, more 
detailed models indeed show that the modified rotation periods give more centrally condensed models for 
Uranus than for Neptune (Nettelmann et al. 2012). 

4 Effect of Giant Impacts 

As mentioned earlier, a striking difference between Uranus and Neptune is Uranus' large axial tilt, and it 
is tempting to speculate that the giant impact which presumably caused the tilt (Safronov 1966) may also 
be responsible for the different internal structures and heat fluxes of the two planets (Stevenson 1986). To 
investigate the effect of the planet's tilt due to giant impacts we adapted a code that follows the motion of 
planetesimals through a protoplanetary envelope (Podolak et al. 1988). Although this code doesn't model 
the details of hydrodynamic effects or the high temperature/density equation of state effects, it can estimate 
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the energy and angular momentum budget when giant impact occurs. For the planetary target we adapt one 
of the models of Podolak et al. (1995) with a rock core surrounded by an ice shell and gas envelope. The 
details of the model aren't important, since we don't know what the pre-impact Uranus and Neptune looked 
like. For our purposes the important thing is that the ice shell extends to 0.75 of the planetary radius which 
is typical of Uranus and Neptune models. At this point there is a sharp density discontinuity between the 
outer envelope and the ice shell. Large planetesimals that penetrate to this region will break up due to the 
sharp increase in ram pressure when they hit the shell. 

The lower panel of Fig. [3] shows the angular momentum deposition by a 5 x 10 3 km rocky body (blue) for 
impact parameters of b = 0.984 Rjjranus (solid) and b = 0.945 Riiranus (dashed), respectively. The angular 
momentum deposition is given in units of Uranus' present angular momentum, while the velocity is given in 
units of Uranus' orbital velocity. The escape velocity from Uranus is roughly three times its orbital velocity, 
so typical encounter velocities should be around iv or ut J$ v < Gv rbit although special circumstances could 
lead to higher/lower velocities. 

As suggested from the figure there are two distinct regimes. For velocities < 3.5v or bit the planetesimal 
has a trajectory that allows it to penetrate into deep enough regions where the ram pressure forces exceed 
the internal strength of the material as well as its gravitational binding, and the planetesimal breaks up in 
the envelope. In that case the planetesimal deposits all of its angular momentum in the planet. For higher 
velocities the planetesimal's trajectory keeps it in the upper envelope where the density is significantly lower 
and the planetesimal simply passes through the planet, causing a sharp drop in the angular momentum 
deposition. As the velocity increases still further the planetesimal leaves more of its angular momentum 
behind as a result of drag forces. For the smaller impact parameter a 5 x 10 3 km planetesimal can deposit 
sufficient angular momentum to change Uranus' tilt without being captured. For a 10 4 km planetesimal (red) 
the behavior is very similar, but the angular momentum deposited in the planet is correspondingly higher, 
and even at an impact velocity of ^v or bit the planetesimal can pass through the planetary envelope and still 
deposit enough angular momentum to account for the tilt. 

The upper panel of Fig. [3] presents the energy (dashed) and angular momentum (solid) deposited vs. 
impact parameter for a 5 x 10 3 km planetesimal impacting the planet. The energy is normalized to the 
gravitational energy difference between a fully mixed planet and a differentiated planet. The angular mo- 
mentum is normalized to Uranus' current angular momentum, and the impact parameter is normalized to 
the planet's radius. The blue and red curves correspond to impact velocities of 5v or bit and 10v or bit, respec- 
tively. It is found that the angular momentum deposited increases as the impact parameter increases until 
some maximum value of b. After that the planetesimal passes through the planet leaving a relatively small 
fraction of its angular momentum behind. The energy deposited at the point of breakup increases as the 
impact parameter decreases, and even for the case of v — 5v or bit is somewhat < 1 in these normalized units 
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so that there's sufficient energy to mix large fractions of the core into the envelope. 

Uranus and Neptune could have both suffered collisions with large bodies. Suppose that Uranus' collision 
was oblique, while Neptune's collision was more radial. In the case of Uranus, the oblique collision could 
not only cause a large tilt in Uranus' spin axis, but could also deposit its energy into a thin shell to inhibit 
convection (Stevenson 1986). A more radial collision for Neptune would mix the interior, and make the 
temperature profile closer to adiabatic, leading to more rapid cooling. This could explain why Neptune has 
a strong internal heat source and Uranus doesn't. It also explains the difference in the cooling histories 
of Uranus and Neptune which was recently found to be even more prominent when using updated model 
atmospheres (Fortney et al. 2011). Indeed, our simulations show that radial collisions, which add relatively 
little angular momentum to the planet can have sufficient energy to mix large fractions of the core, while 
oblique collisions can add large amounts of angular momentum and still allow the core to survive. Thus 
Uranus' large axial tilt may be directly related to the fact that it is more centrally condensed than Neptune. 
This could be tested in the future by more careful measurements of the gravitational moments of Uranus 
and Neptune. 

5 Summary 

The three measured parameters that are used to determine the structure of Uranus and Neptune's (J2, q, u>) 
are not self-consistent and any chosen pair predicts a different inertia factor for the planets. The choice of q 
and J2 which are the ones used in interior modeling, gives very similar inertia factors for Uranus and Neptune 
implying that their internal structures are alike. However, when we use the HAS10 rotation periods, the 
parameters q, f and J2 are more self-consistent and Uranus is found to be more centrally condensed than 
Neptune. As a result, Uranus and Neptune may differ substantially in their internal structures (Nettelmann 
ct al., 2012). 

Two important results derive from our study. First, one must be cautious in drawing conclusions about 
planetary structure simply on the basis of a mass classification. Uranus and Neptune differ by only 10% in 
mass and only 3% in radius, but they may have very different interiors. Second, these differences may be 
caused by secondary effects, such as the impact parameter of a giant impact. These effects have important 
consequences for the classification of intermediate-mass bodies, and more detailed studies are necessary in 
order to evaluate the differences we might expect in the structure and composition of the many intermediate- 
mass cxoplanets. 
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Table 1: Input Parameters 



Parameter 


Jupiter 


Saturn 


Uranus 


Neptune 


Reference Radius 
Period 
/ 

J 2 x 10 6 


71492 
9.925h 
0.06487 
14696.43 ±0.21 


60330 
10.570 h 
0.09796 
16290.71 ±0.27 


26200 
17.240 h 
0.0229 
3341.29 ±0.72 


25225 
16.11h 
0.0171 
3408.43 ± 4.5 



Table 2: Inertia Factor (The letters V and H in parenthesis refer to the Voyager and HAS10 periods, 
respectively.) 



Jupiter 


Saturn 


Uranus (V) 


Neptune (V) 


Uranus (H) 


Neptune (H) 


Method 


0.2661 


0.2278 


0.2527 


0.2067 


0.2293 


0.2549 


q + f 


0.2658 


0.2260 


0.2219 


0.2342 


0.2161 


0.2482 


q + .h 


0.2656 


0.2254 


0.2122 


0.2452 


0.2122 


0.2452 


f + J2 
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Figure 1: Inertia factor vs. rotation period as computed with three different approximations: q and J2 
(maroon), q and / (blue), and / and J2 (green). The upper panel shows Jupiter (solid) and Saturn (dashed). 
The measured rotation periods are given by the vertical black lines for Jupiter (solid) and Saturn (dashed). 
The middle panel shows the same plots for Uranus (solid) and Neptune (dashed) using the Voyager periods. 
The bottom panel is the same, but using HAS10 periods. 
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Figure 2: — J4/J2 vs. normalized core radius for Uranus (red) and Neptune (blue). Solid curves are for 
Voyager rotation periods, dotted curves are for HAS10 periods. The rectangles show the range of models 
with moments that match the observed values. 
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Figure 3: Upper panel: Angular momentum (solid) and energy (dashed) deposited vs. impact parameter 
for v = 5v or bit (blue) and v = lQv or ut (red). Lower Panel: Angular momentum deposited as a function 
of impact velocity for a 10 km radius rock impactor (red) and a 5 x 10 3 km impactor (blue) for impact 
parameters of 0.984 Ruranus (solid) and 0.945 Ruranus (dashed). 



